Simple asymptotic forms for Sommerfeld and Brillouin precursors 
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This article mainly deals with the propagation of step-modulated light pulses in a dense Lorentz- 
medium at distances such that the medium is opaque in a broad spectral region including the carrier 
frequency. The transmitted field is then reduced to the celebrated precursors of Sommerfeld and 
Brillouin, far apart from each other. We obtain simple analytical expressions of the first (Sommer- 
feld) precursor whose shape only depends on the order of the initial discontinuity of the incident 
field and whose amplitude rapidly decreases with this order (rise-time effects). We show that, in 
a strictly speaking asymptotic limit, the second (Brillouin) precursor is entirely determined by the 
frequency-dependence of the medium attenuation and has a Gaussian or Gaussian-derivative shape. 
We point out that this result applies to the precursor directly observed in a Debye medium at 
decimetric wavelengths. When attenuation and group-delay dispersion both contribute to its forma- 
tion, we establish a more general expression of the Brillouin precursor, containing the previous one 
(dominant-attenuation limit) and that obtained by Brillouin (dominant-dispersion limit) as particu- 
lar cases. We finally study the propagation of square or Gaussian pulses and we determine the pulse 
parameters optimizing the Brillouin precursor. Obtained by standard Laplace-Fourier procedures, 
our results are explicit and contrast by their simplicity from those derived by the uniform saddle 
point methods, from which it is very difficult to retrieve our asymptotic forms. 

PACS numbers: 42.25.Bs, 42.50.Md, 41.20.Jb 



I. INTRODUCTION 

More than one century ago, in a short communica- 
tion [l| made at the 79"' congress of the German physi- 
cists, Sommerfeld examined the apparent inconsistency 
between the theory of special relativity and the possibil- 
ity of superluminal group velocity predicted by the clas- 
sical wave theory. Considering an incident wave switched 
on at time t = —T and off at time t = T (square-wave 
nrodulation) , he mathenratically demonstrated that, re- 
gardless of the value of the group velocity at the fre- 
quency of the optical carrier, no signal can be transnrit- 
ted by any linear dispersive-attenuative medium before 
the instant t = —T + z/c, where z is the propagation 
distance and c the velocity of light in vacuum. In the dis- 
cussion following the Sommerfeld's communication, Voigt 
proposed a sinrple physical interpretation of this result. 
He remarked that the front of the wave encounters a 
medium that, due to its inertia, seems optically empty 
and, thus, that the propagation of the very first begin- 
ning of the signal will proceed undisturbed with the ve- 
locity of light in vacuum. In other words, local causality 
implies relativistic causality. The analysis of what hap- 
pens after the arrival of the wavefront was subsequently 
conducted by Sommerfeld and Brillouin in the case of a 
step- wave modulation (field switched on at time i = 0), 
the medium being modeled as an ensemble of damped 
harmonic oscillators with the same resonance frequency 
Wo and the same damping rate 7 (Lorentz medium) [^^ 
|5|. They found that, in suitable conditions, the trans- 
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nritted signal consists in two successive transients (that 
they named "forerunners") preceding the establishment of 
the steady-state field at the frequency Wc of the optical 
carrier (the "main field"). The first and second forerun- 
ners, now called the Sommerfeld and Brillouin precur- 
sors, were associated with the frequencies respectively 
high and low compared to the resonance frequency wq of 
the medium. These results were obtained by means of a 
spectral approach involving the newly developed saddle- 
point method [3j but also classical complex analysis 0] 
and stationary phase method [J|. Following these pio- 
neering works, precursors became a canonical problem 
in electromagnetism and optics [y, 0| • Results complet- 
ing, improving and even correcting those of Sommerfeld 
and Brillouin were obtained by means of uniform asymp- 
totic methods |8l-[ll]. The problem was also studied by 
a purely temporal approach p^l- At the present time, 
the theoretical study of precursors continues to raise a 
considerable interest. An abundant bibliography can be 
found in the recent Oughstun's book \\^. Complemen- 
tary studies on the effects of a finite turn-on time of the 
incident field on the precursors are reported in |14| - [l7| . 

From an experimental point of view, the observation of 
Sommerfeld and Brillouin precursors in the optical range 
raises serious difficulties. Indeed the excitation of the 
Sonrmerfeld and Brillouin precursors requires the cor- 
responding frequencies (respectively high and low com- 
pared to Wo) be present at a significant level in the spec- 
trum of the incident pulse. An experiment intended to 
observe the Brillouin precursor in water is reported in 
|18J . Using pulses at a wavelength of 700 nm with a 
bandwidth of 60 nm, the authors observed pulse breakup 
in a linear regime as well as a sub-exponential attenu- 
ation with distance of the new peak. They attributed 



these features to the formation of a Brihouin precursor. 
This interpretation has been soundly disputed, in partic- 
ular because the pulse bandwidth was in fact not broad 
enough to perform the excitation of precursors [19]. Al- 
ternative explanations of the observations have been pro- 
posed [IQI . 1201 and more recent studies |2ll - l23j have con- 
firmed that a sub-exponential decay of the transmitted 
energy does not prove the formation of precursors. 

While well distinguishable Sommerfeld and Brillouin 
precursors are expected when the medium is opaque in 
a broad spectral region, coherent transients of another 
kind are obtained in the opposite case where the width 
of the opacity region is very small compared to the res- 
onance frequency ujq. They have been naturally named 
resonant precursors |24| but also Sommerfeld-Brillouin 
precursors [25|. Indeed they may be seen as resulting 
from the coalescence of the Sommerfeld and Brillouin pre- 
cursors, originating a well-marked beat when the optical 
thickness of the medium is large enough p6|. The condi- 
tions required to achieve experimental evidence of these 
precursors are relatively easy to meet. They have been 
actually observed in various systems, in particular in a 
molecular gas |2a , in a solid-state sample with a narrow 
exciton line [2^ and in clouds of cold atoms |27l . [28j . 

In the present paper we come back to the study of 
Sommerfeld and Brillouin precursors in a dense Lorentz 
medium, considering the limit where the medium is 
opaque in a spectral region of width large compared to 
the resonance frequency. We remark that these condi- 
tions are met for the parameters considered by Brillouin 
[29] and often referred to in the literature. We then suc- 
ceed in obtaining simple and explicit analytical expres- 
sions of both precursors. When it is necessary, we de- 
termine the range of validity of these analytical solutions 
by comparing them to exact numerical solutions obtained 
by fast Fourier transform (FFT) . The arrangement of our 
paper is as follows. In Section [III we outline the problem 
under consideration and give some general results, useful 
for the following. Section Hill is devoted to the study of 
the Sommerfeld precursor. We establish the correspond- 
ing expression of the impulse response of the medium and 
apply it to obtain a general expression of the precursors 
obtained with causal incident fields. We examine in detail 
the particular cases where the incident field is discontin- 
uous at the initial time or has the canonical form consid- 
ered by Brillouin with eventually a finite rise time. We 
show in Section lTVl that. in a strictly speaking asymptotic 
limit, the impulse response associated with the Brillouin 
precursor is Gaussian and that the Brillouin precursor 
has itself a Gaussian or Gaussian-derivative shape. The 
precursor obtained in a Debye medium is incidentally 
examined. A more general expression of the Brillouin 
precursor in the Lorentz medium is established in Sec- 
tion |Vl containing the previous one and that obtained by 
Brillouin as particular cases. The propagation in both 
media of pulses with a square or Gaussian envelope is fi- 
nally examined in Section IVll and we determine the pulse 
parameters optimizing the Brillouin precursor. We con- 



clude in Section IVIII bv summarizing and discussing our 
main results. 



II. GENERAL ANALYSIS 

We consider a one-dimensional optical wave propagat- 
ing in a Lorentz medium in the z-direction, with an elec- 
tric field linearly polarized in the a;-direction (x, y, z : 
Cartesian coordinates). We denote e(0,i) the algebraic 
amplitude of the field at time t ioT z = (inside the 
medium) and e{z,t) its value after a propagation dis- 
tance z through the medium. The incident field e(0,i) 
being given, the problem is to determine the transmitted 
field e{z,t). We take for e{0,t) the general form 



e(0, t) = u{t) cos(ujct — <^), 



(1) 



including as particular cases the different forms consid- 
ered in the literature, ujc is the frequency of the optical 
carrier, (p is the phase (eventually time-depending) and 
u{t) > is the amplitude modulation or field envelope. 
On the other hand the medium is fully characterized in 
the frequency domain by its transfer function H(z,uj) 
relating the Fourier transform E{z,uj) of e{z,t) to that 
E{0,uj) of e(0,t) [131. 



Eiz,u)^H{z,uj)E{0,uj). 



(2) 



In all the following, we take for t a retarded time equal 
to the real time minus the luminal propagation time z/c 
(retarded-time picture). H{z,uj) then reads 



H{z, Lo) = exp <^ ~i — [n{Lu) — 1] > 



(3) 



Here n{Lj) is the complex refractive index of the medium 
at the frequency w, that is for the Lorentz medium 



1/2 



n{u:) = 1 - ^ 



2i'yuj 



(4) 



where ujq is the resonance frequency, 7 is the damping or 
relaxation rate and ujp is the so-called plasma frequency 
whose square is proportional to the number density of ab- 
sorbers. 3ff [7^(w)] is the usual (real) refractive index n{uj) 
and the absorption coefficient a(w) for the amplitude is 
given by the relation a{u!) = — (w/c)^ [?^(w)]■ 

In the time domain, the medium will be characterized 
by its impulse response h{z,t), inverse Fourier transform 
of H{z,uj), and the transmitted signal e{z,t) is given by 
the convolution product [30] 



e{z,t) = h{z,t)ig)e{0,t). 



(5) 



Some general properties of h{z,t) and e{z,t) can be de- 
duced from Eqs.(I5][S]). First h{z,t) fulfills the condition 
of relativistic causality, namely h{z,t) = for t < |31| . 



Its area reads as /_ h{z,t)dt = H{z,0) = 1. It keeps 



thus constant and normalized to unity regardless of the 
propagation distance z. Consequently i?(z,0) — i?(0,0) 
, that is 



+00 



+00 



e(z, i)dt 



e{0,t)dt. 



(6) 



The area of the optical field (to distinguish from that of 
its envelope) is conserved during the propagation. Fi- 
nally the fact that H{z, 00) = 1 entails that h{z,t) will 
start by a Dirac delta- function S{t). This implies that 
the propagation of the very first beginning of any inci- 
dent signal e(0, t) will always proceed undisturbed at the 
velocity c, in agreement with the Voigt's remark on the 
Sommerfeld's communication [ij. The previous results 
are valid whatever the values of the parameters may be. 
Examine now in what conditions the medium is opaque 
in a broad spectral region. To be definite, we will con- 
sider that the medium is opaque at the frequency tu when 
its optical thickness a{uj)z exceeds 20, the amplitude 
transmission |_ff(z,a;)| = exp [— Q;(a;)z] being then about 
2 X 10^^. Following Sommerfeld [2|, we characterize the 
propagation distance by the parameter ^ = Upz/2c , ho- 
mogeneous to a frequency. For large propagation dis- 
tances 7^/(10wq) ^ 1 and it is easily derived from Eq. 0] 
that the medium will then be opaque in the broad spec- 
tral region w_ < w < w+ with uj+ /tuo ~ ■\/7C/(10i^o) ^'^"^ 
uj-fuo ~ (1+^2/^2)^/^ ^10a;2/(7^). The inequality 
7^ /(lOwg) ^ 1 is over-satisfied for the parameters values 
considered by Brillouin [29], namely tuo — A x 10^^ s^^, 
Wp = 1.24 w^, 7^ = w^/200 and z = lO^^j-^. We then get 
^ = 3.0324 X lO^is-i and 7^/(10a;g) w 5.87 x 10^. Not to 
reduce our study to a particular system or region of the 
spectrum, all the frequencies (the times) will be referred 
in the following to their natural unit luq (I/ujo). 



Figure [T] shows the profiles of the amplitude transmis- 
sion \H(z,uj)\) = exp [— a (w) z] as a function of the re- 
duced frequency lo/luq in the Brillouin conditions (curve 
a) and for propagation distances 10, 100, 1 000 and 10 
000 times shorter (curves b to e). 

The medium being opaque for aj_ < a; < aj_|_, the trans- 
fer function may be written as 



H{z,uj) = Hs{z,uj) + Hb{z,uj), 



(7) 



with Hs{z,uj) « for lo < aj_|_ and Hb{z,uj) w for 
uj > aj_. Hs and Hb are respectively associated with 
the Sommerfeld and the Brillouin precursor. For cj = 0, 
Hs{z,0) w and Hb{z,Q) w H{z,{)) = 1. As long as Uc 
lies in the opacity region, this implies that the Sommer- 
feld precursor will have a zero area while the area of the 
Brillouin precursor will he equal to that of the incident 
field. 

The formation of the optical precursors is generally 
governed by combined effects of attenuation (consid- 
ered above) and dispersion. The dispersion effects can 
be soundly characterized by the group delay Tg{z,uj) — 
—d^/duj = z/vg{oj) — z/c, where $(2:, uj) is the argument 
of H{z, ijj) and Vg{u)) the group velocity [31j. We remark 
that the regions of anomalous dispersion {dn/diu < 0) 
or of superluminal group velocity (r^ < ) has a width 
smaller than luq and are entirely comprised inside the 
opacity region. The corresponding frequencies will thus 
not directly contribute to the formation of precursors. 
For the high and low frequencies respectively associated 
with the Sommerfeld and Brillouin precursors, we get the 
asymptotic forms Tg 
where 



e/w2 P and Tg « tB+Loy{r,b^) 



tB = 



[n{0) - 1] z _ 2^ 



OJ^ 




(8) 




10"' 10"' 1 10' 10' 
Frequency (in cOg units) 

Figure 1: Amplitude transmission \H{z,lj\ of the medium as 
a function of tlie frequency modulus \uj\ (logarithmic scale). 
Parameters (in luq units) : LUp = 1.11 , 7 — 0.0707 and ^ — 
8.31 x 10^ for the curve (a) corresponding to the Brillouin 
choice (z — 10~^m ). The curves (b), (c),(d) and (e) are 
obtained for propagation distances (and thus ^ ) respectively 
10, 100, 1 000 and 10 000 times smaller. 
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1/6 



/ 1 + 



(9) 



(10) 



tB — Tg(z,0) — Tg(z, 00) is obviously indicative of the 
time delay of the Brillouin precursor (low frequency) with 
respect to the Sommerfeld precursor (high frequency). 
The two precursors will be fully separated when ts is 
much larger than the damping time I/7. Since jts — 
O (71^/wo) 7 this condition is automatically fulfilled when 
the condition of broad opacity- region [7^/(10a;o) ^ l] 
holds. Another important point is that Tg is minimum 
(stationary) for uj ^ 00 and w — >■ 0. As pointed out by 
Brillouin [4] , this ensures that the precursors will not be 
washed out by the group velocity dispersion. 



III. SOMMERFELD PRECURSOR 



product of Eq.® takes the form: 



A. Transfer function Hs{z,u)) and impulse response 



In the hmit considered here w^ > w^ 3> Wg and 
Hs{z, w) takes the following asymptotic form, accounting 
for both dispersion (main contribution) and attenuation. 



Hs{z,u) w exp 



^ 



27 



(11) 



The corresponding impulse response hs{z,t) is easily de- 
termined by using standard results of Laplace transforms 
[Hi. We get 



hsiz,t)^d{t) 



#.( 



W^t 



-2'yt 



UH(t), (12) 



where Jn(s) and VLnit) respectively designate the first 
kind Bessel-function of index n and the Heaviside unit- 
step function. Except for their very first oscillation, 
the Bessel functions J„(s) are perfectly approximated by 
their asymptotic form 



Jn(s) 



COS I 



s — n 

2 4 



(13) 



es{z,t)= / hs{z,e)e{Q,t~e)de, (15) 



that can be transformed by repeated integrations per 
parts to yield 



es{z,t) 



E 

n=0 



drM^+'\z,t). 



(16) 



Here dn is the discontinuity of the n**^ derivative of e(0, t) 
at the initial time |33|| and /^"■'(t) is a short-hand nota- 



tion for /__^ /_'^ • • • /_'^' f{tn)dtn ■ ■ ■ dt2dti. In a fre- 



quency description, the previous result can be retrieved 
by expanding the Fourier transform E{0,lu) of e(0,i) in 
powers of l/ito and exploiting the equivalence between 
multiplication by l/icu in the frequency domain and inte- 
gration in the time domain [30J. Writing the impulse re- 
sponse under the form hs{z,t) = ks{z,t)exp{—2jt), we 
easily show by means of standard Laplace procedures |32| 

that 4"^^^(z,i) = (i/0"/^Jn(2^/Ct)uH(i). Insofar as 
ks{z,t) is very rapidly varying compared to exp(— 27t), 

hg {z,t) ~ kg {z,t) exp(— 27i) and we finally get 



and the impulse response hs{z,t) can be characterized 
by an instantaneous frequency oj ^ d (2^/0) /dt = 

\/£,/t. The range of validity of Eq. (fT2l) may be esti- 
mated by determining the change SHs{z,ui) of Hs{z,u)) 
due to the first term neglected in the asymptotic ex- 
pansion of \n\Hs{z,ijo)] used to obtain Ea. (|ll|) . We 
find 6Hs{z,Lu)/Hs{z,uj) — O (^Wq/cj^), negligible when 
w^ > ^uj^^ i.e. when ^^/^ > uj^t^^^. In fact, Eg. p^ 
fits very well the exact impulse response as soon as ^^/^ 
exceeds Wgi"^/^ by a factor a/IO (half an order of magni- 
tude). This is achieved as long as i < ^5, with 



^0*5 = 



^ 



lOwo 



(14) 



In a strict asymptotic limit (z — ^ cxi), ^5 — >■ cxd and 
exp (—27^5) — > 0. As expected, the entirety of the im- 
pulse response is then reproduced by Ea. (|T2|) . 



B. Precursor originated by a causal incident field 

The Sommerfeld precursor es{z,t) is obtained by 
convoluting hs{z,t) with the incident field e(0,i) = 
u{t) cos {uJct — if) introduced in the general analysis 
[Eq.([T|)]. We are mainly interested here in the physi- 
cal case where the incident field is causal [e(0,i) = for 
t < 0], u{t) being either a unit step u//(i) or a function 
monotonously rising from to 1 with a rate r < Uc for 
t > (step or step- like modulation). The convolution 



es 



n=0 



ri/2 



(z,t)«^d„(-j J„(2Ve^)cxp(-27i)uff(i). 



(17) 

The n*'' term of the series has a maximal amplitude ao = 
Idol at t = to = for n = and 



1 



Mnl 



2r7.-l 
8e 



(2n-l)/4 



1/4 



ilO 



-n/2 



(18) 



at t w tn = (2?! - 1) 787 for n > 0. Since ^ oc z, Eq.lHl 
shows that, for large propagation distance, a„ rapidly de- 
creases with n, so that a good approximation of the exact 
result is obtained by keeping only the first term n — p 
of the series for which dp y^ 0. In the frequency descrip- 
tion, this amounts to restrict the asymptotic expansion 
of E{0, oj) to its first non vanishing term (7|. We then get 



es{z,t) ^dp[- 



p/2 



Jp(2Ve0exp(-27t)uH(i). (19) 



Denoting q is the next integer following p for which 



dg 7^ 0, Eq. (jl9p is exact when e — Oq/op w 0, and 
exp (—27^5) sa 0. These conditions are met in the strict 
asymptotic limit and closely approached for the propaga- 
tion distance considered by Brillouin. At distances that 
may be 1 000 times smaller (simple asymptotic limit); 
we shall see that Eq. p^ enables us to correctly repro- 
duce the essential features of the precursor originated by 
representative incident fields. 



Precursor originated by a discontinuous 
incident field 



We consider first the instructive case where e{0,t) = 
Unit) cos (Wfji) for which p = with do = 1 lS3 and q = 2 



and q ~ 3 with d^ 



We then get 



with ^2 



Eg. p^ then reads as 



es(z, t) w Jo(2 VCO exp (-27^) u^ (t), 



(20) 



with £ « 0.13a;27-3/4^-5/4 [g^g gq_ ^ r^^ie precursor 
does not depend on ujc and the initial discontinuity of 
the incident field is integrally transmitted, in agreement 
with the general analysis. For Wc < uj+ = vtC/IO (opac- 
ity condition), e is always smaller than 0.013 (7/O , 
that is about 2.2 x 10~* in the Brillouin conditions and 
1.2 X 10~^ for a propagation distance 1 000 times smaller 
(simple asymptotic limit). In the first case, uj^ts = 44 
and exp (—27^5) « 2 x lO^'^. As previously indicated, 
we are then close to the strict asymptotic limit and the 
precursor is perfectly reproduced by its asymptotic form 
at any time where it has a significant amplitude. 




5 10 15 

Retarded Time (in I/coq units) 

Figure 2: Sommerfeld precursor originated by the incident 
field cos(a;ct)u£r(t). The solid (dashed) line is the exact nu- 
merical solution (the approximate analytic solution). Param- 
eters (in Wo units) : u)c — 1, u)p — 1.11, 7 — 0.0707 and 
^ = 831. Inset : enlargement of the tail of the precursor. 



This remark also holds for the cases considered in the 
following subsections. In the simple asymptotic limit 
LijQts — 4.4 and, as expected, Eq. (|20p perfectly fits the 
exact solution for Lo^t < 4.4. For larger times, the fit 
remains very good except for a slight drift of the instan- 
taneous frequency of the oscillations whose envelope is 
very well reproduced at any time (Fig|2]). 



es(z,i) fnuj, 



:j|ji(2v/^)exp(-27i)uH(i), 



(21) 



with e w 0.34(a;^/7^). The result given Eq.(l2l]) dif- 
fers from that originally obtained by Sommerfeld [2] by 
the presence of the damping term exp (—27^) . Though 
the formation of the Sommerfeld precursor is mainly gov- 
erned by the medium dispersion, the presence of this term 
(associated with the absorption) is obviously necessary 
to avoid that esiz^t) diverges with time. The precursor 
attains its maximum at i « ii = 1/(87) ("^o^i = 1-77 
) and its amplitude as = ai ~ 0.26 lOc^^^^^^^'^^'^ is 
proportional to ujc- For lOc = wq, as ~ 1.8 x 10^^ 
with e « 5.8 x 10^^ in the Brillouin conditions whereas 
as ~ 3.25 X 10"^ with e « 5.8 x 10~^ in the simple 
asymptotic limit. In the latter case. Fig [3] shows that 
Eq. (PT|) actually fits very well the exact result for t < ts, 
again with a slight drift of the instantaneous frequency 
of the oscillations for t > ts- In order to check the pro- 
portionality of the precursor to Wc, we have compared 
the exact forms of (uo/ujc) Gs{z,t) obtained when ujc 
lies at the boundaries w_ or lo^ of the opacity region 
to that obtained when Wc = wo- As expected we have 
found that the three results are nearly undistinguish- 
able, except for an amplitude 1.3% larger for lOc — ^+ 
(below the corresponding value of e, namely e — 0.034). 
For this value of Wc, the amplitude of the precursor is 
as w 0.082 (7/0^^'*, that is 1.40 x 10"^ in the Brillouin 
conditions and 7.9 x lO^'' in the simple asymptotic limit. 



xlO^ 



Retarded Time (in I/cOq units) 

Figure 3: Sommerfeld precursor originated by the canonical 
incident field sm{ujct)uH{t)- The solid (dashed) line is the 
exact numerical solution (the approximate analytic solution). 
Parameters as in Fig[21 



D. Precursor originated by the canonical incident 
field of Sommerfeld and Brillouin 



Following Sommerfeld and Brillouin, most authors 
have considered an incident field of the canonical form 
e(0,t) = u nit) sin (ujct) for which p = 1 with di — ujc 



E. Rise-time effects 

A gradual turning on of the incident field is ex- 
pected to reduce the amplitude of the Sommerfeld pre- 
cursor. To study this so-called rise-time effect, Cia- 
rkowski [ij,[i3| has considered the incident field e(0, i) = 



tanh (ri) sin (ujct) unit) whose envelope has a 10 — 90% 
rise time T^ ~ 1.37/r. In this case p — 2 with ^2 = 2ra;c, 



q = 4 with d4 = -4wcr (2r"^ -f 
form of the precursor reads as 



and the asymptotic 



simulations for e(0,i) = sin (wci) [1 + erf (ri)] /2 where 
erf(s) designates the error function. For z and r = loq 
as in Fig. 4, we get p w 1.4 x lO'' instead of p w 7.7 for 
e(0,t) — tanh(ri) sin{uJct) unit). 



es{z,t) « 2u)cr 



-)j2(2v1i)exp(-27t)uff(t), (22) 



with £ « 1.21 (2r2 



i^c) /tC- The precursor attains 
its maximum at i w ^2 = 8/(87) (^^0^2 ~ 5.3) with 
an amplitude as — 0-2 ~ 0.26 rtOc^^'^^^^,^^^^- Com- 
pared to the precursor obtained with the canonical in- 
cident field [Eg. ([21]) ]. the maximum is shifted to larger 
time (i2 = 3ii) and its amplitude is reduced by a fac- 
tor p sa VT^/r- Fig|31 obtained in the simple asymptotic 
limit, shows that Eg. ((22)) fits quite satisfactorily the exact 
precursor though its maximum now lies at a time slightly 
larger than 1$- To check that the precursor is mainly de- 
termined by the lowest order initial discontinuity of the 
incident field regardless of its subsequent evolution, we 
have compared the precursor obtained when the enve- 
lope tanh (rt) uuit) is replaced by (1 — e"""*) unit), hav- 
ing the same initial discontinuity. Though q = 3 (instead 
of 4) and Tr ~ 2.20/r (instead of 1.37/r), we have found 
that the precursor is actually very close to the previous 
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Figure 4: Sommerfeld precursor originated by the incident 
field e(0, t) = tanh(rt) sin(ajci) unit)- The solid (dashed) 
line is the exact numerical solution (the approximate analytic 
solution) obtained for r = ojq. Other parameters as in Fig[2] 

Other things being equal, the reduction of the ampli- 
tude of the precursor is more and more important when 
the incident field is applied more and more smoothly, 
that is when the order p of its initial discontinuity in- 
creases. It is easily deduced from Eq. (|18P that for p > 2, 

p = [(V7?/r)^"^] oc (T^v^)^"\ At the light of this 

result, dramatic rise time effects are expected when the 
incident field is ideally smooth, i.e. analytic with con- 
tinuous derivatives in every point. Such fields have been 
considered [ij, [la, |33[ though they are not causal and, 
strictly speaking, not physically realizable (in the sense 
of the linear systems theory) . We have made numerical 



IV. BRILLOUIN PRECURSOR IN THE STRICT 
ASYMPTOTIC LIMIT 

A. Transfer function Hb{z,(jj) and impulse response 



In the limit considered now ui'^ < (jPl ^ Wq 
Hb{z,lj) is conveniently developed under the form 



and 



Hb(z,lo) — exp 



E 



{-iujf' 



■kn{z) 



(23) 



Here kn{z) are the so-called cumulants, generally in- 
troduced in probability theory |32[ | , b ut also quite use- 
ful to study deterministic signals |35l [36| . The cumu- 
lants ki{z), k2{z) and ^3(2) have remarkable proper- 

ties, fci(z) and fcj (^) respectively are the center-of- 
mass and the root-mean-square duration of the impulse 
response hsiz, t), inverse Fourier transform of Hb{z,u}), 
whereas k{z) = k^{z)/kr, (z) is its normalized asym- 
metry or skewness [32|. From Eqs. ()3l4p . we easily get 
ki = Ib (as expected), /c2 = 47/(35^) , /C3 = -2/(776^) 
and K — —(1/477) (3^/7) , where t^, h and 77 are defined 
by Eas.([5] fTn[) . When z — >■ 00 (strict asymptotic limit), 
K ex b'^/'^ ex z^^/"^ — >■ and the expansion of Ea. ([23[) may 
be limited to the term n = 2. Taking a new origin of 
time at i = is, the transfer function then reads as 



(24) 



i/B(z,cj)«exp(- — 



where /3 = y'SUySj oc 1/%/z is very small compared to 
Wo- This Gaussian form is that of the normal distribution 
derived by means of the central limit theorem in prob- 
ability theory. This theorem can also be used to obtain 
an approximate evaluation of the convolution of n deter- 
ministic functions [30|- It can be applied to our case by 
splitting the medium into n cascaded sections, hB{z,t) 
being the convolution of the impulses responses of each 
section. By calculating the inverse Fourier transform of 
Hb{z,uj), we get 



hB{z,t) 



—j= exp (- 



-P^t 



") 



(25) 



where t' = t — Ib- The impulse response has a duration 
(amplitude) proportional (inversely proportional) to y^, 
with an area constantly equal to 1 (in agreement with 
the general analysis). We remark that the approxima- 
tion leading to Ea. ([^^ and Ea. ([^5[) . valid in the strict 
asymptotic limit, amounts to neglect the effects of the 
group delay dispersion, the formation of the Brillouin 



precursor being then governed by the frequency depen- 
dence of the medium attenuation (dominant-attenuation 
hmit) . 

The Gaussian forms of Eq.(|24|) and Eq. (|25l) are not 
specific to the Lorentz medium but have some generahty 
|37J . They hold for the Debye medium [33 , for some ran- 
dom media [39] and, more generally, whenever the trans- 
fer function of the medium can be expanded in cumulants 
and the propagation distance is such that |k| <^ 1. Stoudt 
et al. [38] showed in particular that the results of their ex- 
periments on water (Debye medium) at decimetric wave- 
lengths can be numerically reproduced by neglecting the 
group delay dispersion, as it has been made to obtain 
Eq. (l24l) . See also |4(]| - [43 |. Using a purely temporal ap- 
proach, Karlsson and Ritke [12| early remarked that the 
impulse response of the Debye medium is very close to a 
normalized Gaussian. This property is obviously a conse- 
quence of the previous analysis. The complex refractive 

index now reads as n{u}) = [l + («o ^ l) / (1 + ^i^t)] 
where no is the refractive index for cj — > and r 
is the relaxation time for the orientation of the polar 
molecules [4lj. Including n(w) in Eq.([3]) and following 
the procedure used for the Lorentz medium, we easily 

1/9 

get (3 = [2 (tIq — l) rz/crio] and, taking into account 

^ 1, K « 2.25-\/cr/noz. Note that /3 and n de- 



is the algebraic area A of the incident field (see Sec. Illl) . 
we finally get: 



that 



pends on z as l/^/z (as in the Lorentz medium). The nor- 
malized Gaussian of Eq. ipS)) will thus also be obtained for 
sufhcie nt p ropagation distances. Using the parameters of 
water |4l|, namely ng = V79 and t — 8.5 x 10~-^^s, we 
find that the skewness of 5.2%, obtained in a Lorentz 
medium for a propagation distance larger by more of 
four orders of magnitude than the optical wavelengths 
considered, is now attained for a propagation distance 
z « 0.55 m comparable to the wavelengths involved in 
the experiments reported in [38;| . Despite strongly differ- 
ent scales, Brillouin precursors in the Lorentz medium in 
the strict asymptotic limit and in the Debye medium per- 
tain to the same physics, namely that of the dominant- 
attenuation limit, and will be described by the same laws. 
On the other hand, the Debye medium is fully opaque 
at high frequency and Sommerfeld precursors cannot be 
generated in this medium. 



B. Precursor generated by an incident field of 
non-zero area 

The Brillouin precursor generated by an arbitrary in- 
cident field e(0,t) is obtained by convoluting the lat- 
ter with /iB(z,i) or by multiplying its Fourier trans- 
form E{0,uj) by Hb{z,uj) and determining the inverse 
Fourier transform of the product. We consider first the 
case where e(0, t) is rapidly varying compared to hsiz, t). 
This requires in particular that uJc ^ P- Compared 
to _E(0,aj), Hb{z,uj) then appears as a narrow peak 
centered on w = and, provided that £'(0,0) 7^ 0, 
Eb{z,w) w £(0,0) Hb{z,uj). Remembering that £(0,0) 



es 



(z, t) « AhB{z, i) = 4= «^^P (-/3^i'^) • (26) 



For the canonical incident field sin (w^i) u/f (t), £(0, 0) = 
1/ljc and the precursor has an amplitude as = 
/3/ {uJc\/t^) inversely proportional to ujc (no matter its 
value provided that Wc ^ /3) and to y^z. Note that the 
law aB oc 1/v^) sometimes considered as general, is only 
valid in the strict asymptotic limit considered here (for 
which \k\ <C 1). Fig[S]shows that the precursor obtained 



0.01 




664000 



665000 



666000 



667000 



668000 



Retarded Time (in l/m. units ) 



Figure 5: Brillouin precursor obtained in the Brillouin con- 
ditions, namely for ujc = 0.1, cjp = 1.11, 7 = 0.0707 
and ^ = 8.31 x 10^ (in ujq units). For these parameters, 
ujotB ~ 6.654 X 10^ and /3 ^ 1.78 x lO'^wo = 1.78 x lO'^Wc. 
Solid lines (bullets •) are the exact numerical solutions (the 
analytic solutions). Curve (a) is the precursor obtained 
with the canonical incident field sai{u}ct)uH(t). The precur- 
sor of curve (b) is originated by the incident field e(0, t) = 
sin(ajci) [1 + erf(ri)] /2 for r — ljc/2-\/2 . Inset: Sommerfeld 
precursor obtained in the conditions of curve (a). It fully 
vanishes in the conditions of curve (b). 

in all the Brillouin conditions [curve (a)] is perfectly fit- 
ted by the Gaussian form of Eg. ip^ . We incidentally 
note that, for the carrier frequency retained by Brillouin 
{lOc = Wo/10), the medium is fully opaque at this fre- 
quency \a (wc) z sa 800], in contradiction with his artist's 
view showing a "main field" (at Wc) larger than the pre- 
cursors. On the other hand, the condition Wc ^ /? is 
well satisfied. The inset in FigjS] shows the Sommerfeld 
precursor obtained in the same conditions. As already 
mentioned, it is perfectly fitted by the analytical expres- 
sion of Eq. (|2ip . Note however that its amplitude is about 
four orders of magnitude smaller than that of the Bril- 
louin precursor. Eg. ipS)) also holds when the envelope 
of the incident field rises in a finite time provided that 
the rate r , as Wc, is large compared to j3. Curve (b) 
of Figini shows the Brillouin precursor generated by the 
incident field e(0, t) — sin {wct) [1 -I- erf {rt)] /2. We have 
then £(0, 0) — (I/^c) exp (— w^/4r^) and the area of the 
incident pulse, equal to \/uJc for r — >■ cx), falls to l/2a;c 
for r = ijJcl2\J\i\ (2) (r w O.GOwc)- As expected, the 



Brillouin precursor is identical to the previous one with 
ampHtude reduced by half and the corresponding Som- 
merfeld precursor completely vanishes. 



Precursor originated by an incident field of zero 
area 
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Figure 6: Brillouin precursor obtained witii the incident fields 
(1 — e"""*) cos(ajcf)uH(t) for (a) r — >■ oo, (b) r = GSoJc and 
(c) r = 20a;c (solid lines). Other parameters as in Fig[5] 
The bullets • correspond to the analytical solutions given by 
Ea. (|27|l or by the combination of this equation with Ea. (|26p . 



Even if Wc,?' ^ /?, Eq. ((26|) obviously fails when A = 
E{0, 0) — 0. This occurs in particular in the extreme 
case where the incident field is instantaneous turned on, 
with e(0,t) = cos (wci) u/f (t) . It is then necessary to 
consider the next term in the expansion of E(0,uj) in 
powers of iu. We get in this case E{Q,uj) « ioj/iJl and 
Eb{z,uj) ~ iijjHB{z,Lo)/uj^. Using the correspondence 
iuj d/dt between frequency and time descriptions [30] 
and denoting by a dot the time derivative, we finally get: 



eB{z,t) w -—hB{z,ujt) 



2/32 



/3i'exp(-/32i'2) 



(27) 
As shown FiglH] [curve (a)], the analytical expression of 
Eq. (|27p perfectly fits the exact numerical results obtained 
by FFT. The precursor is a Gaussian derivative with a 
peak amplitude ub — [2/(7re)] {fi/uic) , smaller than 
that attained with the canonical incident field by a fac- 
tor ajcA/e/(/3-\/2) (~ 65 in all the Brillouin conditions) 
and decreasing much more rapidly with the propaga- 
tion distance (as 1/z instead of as Xj^/z). We how- 
ever remark that the precursor so obtained is not ro- 
bust. Indeed it suffices that the incident field suffers a 
short rise time to retrieve a precursor mainly governed 
by the area law of Ea. ((26|) . To illustrate this point, 
we have again considered an incident field of the form 
(1 — e"''*) cos (wci) u/f (t) that tends to cos (wct) Uij(i) for 
r — >■ oo. For r '^ lOc (very short rise time), £'(0,0;) k, 



— 1/r -I- iujJLij^. The incident field has gained a (nega- 
tive) area A — ^1/r. The precursor is then the sum 
of two contributions, respectively given by Eq. (j26P with 
A = -I/r and by Eq.dl?]). Curve (b) of FigE] shows 
the result obtained when the two contributions have the 
same amplitude, that is when r /ujc ~ '^c\/e/(/?v2) ~ 65 
. When r decreases by remaining large compared to Wc, 
the Gaussian part of the precursor rapidly prevails on 
the Gaussian-derivative part and, as shows [curve (c)], 
the precursor becomes nearly Gaussian (downwards) for 
r as large as 20a;c. 



D. Case where the carrier frequency lies below the 
opacity region 

The previous results are valid for the Lorentz medium 
in the strict asymptotic limit (also as in the Debye 
medium) when Wc 3> /?, that is when lOc lies in the opacity 
region. Fortunately enough, the simplicity of the Gaus- 
sian impulse response enables us to obtain exact expres- 
sions of the transmitted field for arbitrary values of the 
ratio uJcl P- This occurs in the Lorentz medium when Wc 
resides below the opacity region and direct observations 
of the field transmitted in such conditions have been per- 
formed by Stoudt et al. in a Debye medium [38). The 
transmitted field e{z,t) is calculated directly in the time 
domain by convoluting hB{z,t) given Eq. (l25[) with the 
incident field. For the canonical incident field, the con- 
volution product can be written as: 



e(z,i) 



/3 



-'^'^' sin [coc [f - 9)] d9. 



(28) 



After some simple transformations, we finally get 



e{z,t) = ie-"'/4^'3 



1 + erf 



^2/3 



Aujct' 



e 



(29) 
where e~'^<:/'^^ « Q-a(uja)z g^j^^^ g^g previously, t' = t — tB- 
For t' -^ oo, e(z,i) tends to e~'^o/^P sin (wci') which is 
nothing but that the steady state or main field, not neg- 
ligible when ujc and j3 are comparable. If we take Ib 
(I//3) as time origin (time unit), the transmitted field 
only depends on the ratio uJc/ P, regardless of the par- 
ticular system considered. When Uc ^ /3, it tends to 
/3/ (ajc\/7r) exp (— /32<'2) in agreement with Eq. ipSI) . the 
main field being then negligible. When ujc > 4/3, Eq. ([29|) 
is well approximated by the expression: 



e(.,i)«i±^^sin(c.,i')e-"(-)^ 



LOcX^Tl 



(30) 



where /3' = /3 (l-f 2(3^ /ul) -^ ji for w^ > P- The first 
(second) term of Eg. ipH)) obviously corresponds to the 



main field (the Brillouin precursor). Figure [7] shows the 
transmitted field as a function of /3i' — /3{t — tB) for 
ujc ~ 3.84/3 and ujc ~ 7.67/3 (inset). In the study on wa- 
ter (Debye medium) at decimetric wavelengths [38], these 
values are obtained with utr — 2ttx lO^s^^, for z — 0.75m 



and z = 3 m respectively. As expected Eq. ([291) perfectly 
fits the exact numerical result in both cases. Eq. (|30p pro- 
vides a good approximation for Wc ~ 3.84/3, excellent for 
uJc ~ 7.67/3. In the latter case, the Brillouin precursor 
prevails over the main field whose relative amplitude is 
negligible. The signals shown FiglT] are in good agree- 
ment with those directly observed in the experiments re- 
ported in j38i] . 
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Figure 7: Brillouin precursor and main field obtained for 
UJc ~ 3.84/3 as a function of /3(t — ts). The solid line, the 
bullets • and the dashed line are respectively the exact nu- 
merical solution, the analytical solution given Ea. p9|l and its 
approximate form given Ea. (|30p . Inset: Brillouin precursor 
obtained for cjc ~ 7.67/3. The two analytical solutions are 
undistinguishable in this case and the amplitude of the main 
field is negligible. 



EXTENDED EXPRESSION OF THE 
BRILLOUIN PRECURSOR 



We come back in this section to the Brillouin pre- 
cursor in the Lorentz medium. Numerical simulations 
show that the solutions obtained in the strict asymptotic 
or dominant-attenuation limit continue to provide good 
(not too bad) approximations of the exact solutions when 
the propagation distance is 10 times (100 times) shorter 
than that considered by Brillouin [29|, though the skew- 
ness K then rises up to 16% (52%). For shorter distances, 
it is obviously necessary to take into account the effects of 
the group-delay dispersion neglected in the strict asymp- 
totic approximation. 



A. Transfer function Hb{z,uj) and impulse response 



Taking into account the term in w^ in Eq. ([231 
transfer function then reads as 



the 



Hb{z,uj) sa exp 



-iujIb 



Zrjb^ 



(w^ - 2irijuj'^) 



, (31) 



where ts, b and 77 are defined by Eas.([5] fTU)) . with 



27/(353) = 1/4^2 Remarking that 



2ir]juj ] is the 



beginning of (lu — 2ir]j/3) and taking a new origin of 
time at ts + 4777^/96^, we get: 



Hb{z,uj) « exp 



3,63 V""^'^^ 



772/27 
3 V36 



(32) 
By means of an inverse Fourier transform, we finally find: 

/iB(z,i) K BAif-^rj^^^btA exp(-2777i"/3). (33) 



Here B 



V 



1/35 



exp 



- (77 V3) (27/36)^ 



f' 



t-tB 



47772/953 and Ai(s) designates the Airy function. The 
range of validity of Eg. ([55)1 can be roughly estimated 
by means of a strategy similar to that used for the 
Sommerfeld precursor. By taking account of the cu- 
mulants k^ (correction of the attenuation) and /cs (cor- 
rection of the dispersion), the transfer function associ- 
ated with the Brillouin precursor approximately reads as 
Hb{z,uj) X (1 — a4UJ^ — ia^uj^^ where 04 = — A:4/24 > 
and 05 — fc5/120 > . Hb{z,uj) will be a good ap- 
proximation if a4(u;^ and 050;^ are small compared to 
1 (say < I/a/IO). For sake of simplicity, we take for 
the ratios ujp/ujQ and ^ /loq the values retained by Bril- 
louin, representative of a dense Lorentz medium with 
moderate damping. We get then 77 « 1.018 « 1. Be- 
sides, in a cavalier manner, we assimilate uj to the in- 
stantaneous frequency derived from the asymptotic form 



Ai(- 



-l/2„-l/4 



s-i/"* sin (2s3/2/3 -|_ 7r/4) that provides a 



good approximation of Ai(— s) when s > 1. We get so 
\/b^t" . With all these hypotheses, we finally find 



UJ 



that the corrections due to the cumulants k^ and k^ will 
be small ii ujQt" < 2{ujo/b) and ujQt" < {ujo/b) , re- 
spectively. Despite the roughness of the procedure lead- 
ing to these conditions, it will appear below that they are 
realistic and even too severe. 



B. Precursor generated by the canonical incident 
field 



When hB{z,t) is slowly varying compared to e(0,i), 
the Brillouin originated by the canonical incident field 
sin (Wfji) u// (t) takes again the simple form 6^(2, t) — 
A hB{z, t), that is 



eB[z,t) w — At 

UJc 



(-^Tj^^^btA exp (-2777^/3) . (34) 



10 



It is assumed by writing Eq. p4p that the instantaneous 
frequency vWW is small compared to ujc (say VWW < 
Wc/a/IO) and that the conditions of validity of hB{z,t) 
are met. All these restrictions are summarized by the 
inequality 



UJQt" < niin 



2 {cjo/bf^ , {uJo/bf , ^o^.VlOfo'l . (35) 



Fig[5] shows the Brillouin precursor obtained in the sim- 
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Figure 8: Brillouin precursor obtained in the simple asymp- 
totic limit with the canonical incident field sin(ajct)'"£r(t). Pa- 
rameters (in ujo units): cjc = 1, cJp = 1.11, 7 = 0.0707 and 
^ = 831, leading to ujots ~ 665.4, 6 ^ 8.44 x 10"^a;o and 
P « 5.64 X IO'^Wq. The solid line, the bullets • and the 
dashed line are respectively the exact numerical solution, the 
analytical solution given Eq. (|34|) and the Gaussian that would 
be obtained in the dominant-attenuation approximation. The 
conditions are those of Fig[3l The corresponding Sommerfeld 
precursor is given in inset for reference. 

pie asymptotic limit considered in the study of the Som- 
merfeld precursor (FiglS]). The inequality of Ea. (|35p then 
leads to tuot < min [750, 760, 840] . Insofar as the am- 
plitude of the precursor is negligible for ujQt = 750, the 
analytical expression of Eg. dMj) perfectly fits the exact 
numerical result. 

Surprisingly enough, Eq. (j34p remains a not too bad 
approximation of the exact result even when the opacity 
region is not broad in the sense given to this expression 
in the present paper. Fig |9] shows the precursor obtained 
at a distance ten times smaller than the previous one. 
Though the width of the opacity region is then of the 
order of ojq [see curve (e) of FiglT], the entirety of the 
first oscillation of the Brillouin precursor is very well fit- 
ted by Ea. p4p . The corresponding Sommerfeld precur- 
sor (inset) is itself well reproduced by Ea. ([?T|) up to its 
maximum. 



C. Dominant-dispersion limit 

The expression of the Brillouin precursor given by 
Eg. (1341) obviously includes as particular case the Gaus- 
sian obtained in the dominant- attenuation limit. In fact, 
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Figure 9: Comparison of the Brillouin precursor obtained out- 
side the asymptotic limit (solid line) with the analytical forms 
given Eq. (|34|) (•) and Ea. (|36p (dashed line). Parameters (in 
LJo units): uuc = 1, tOp = 1.11, 7 = 0.0707 and ^ = 83.1, lead- 
ing to LJotB ~ 66.54, b ~ 0.182a;o and /3 ~ 0.178a;o. Inset: 
corresponding Sommerfeld precursor (solid line) compared to 
the analytic form given Eq.(|21[I (dashed line). 



retrieving the Gaussian precursor directly from Eq. ([M|) 
requires long and tedious calculations and this probably 
explains why the Gaussian solution has been generally 
overlooked. An other particular form of Ea. (p4)) . also of 
special importance, is that obtained when the damping 
is very small, so that the formation of the Brillouin pre- 
cursor is mainly governed by the group delay dispersion 
(dominant-dispersion limit). This requires in particular 
that 7 < 6. We then get <" « t - ts, S « 6 and 



65(2, t) w — Ai [—b {t — ts)] exp 



2 , 

^37(^ 



ts) 



(36) 



Except for the exponential damping term, this result was 
established by Brillouin himself by means of the method 
of stationary phase |J, |4J[ . When the group-delay dis- 
persion is fully dominant (say when 7/6 < l/lOO), the 
precursor has a well nrarked oscillatory behavior with 
a very weak damping and its maximum practically co- 
incides with the first maximum of Ai [— 6 (i — is)], at- 
tained for i — ts w l,02/&. The corresponding ampli- 
tude is Ob ~ 0.536 (6/cjc) that scales as z"^/^, instead 
of as z~^'^ in the strict or dominant-attenuation limit. 
FiglTU] shows an example of Brillouin precursor obtained 
in such conditions (7/6 w 3.9 x 10~^). It is worth em- 
phasizing that, since h a z^^''^, the condition 7/6 ^ I 
requires that the propagation distance is not too large. 
On the other hand, it should be large enough for the in- 
equality of Eg. ([55]) to be satished for a time larger or 
at least comparable to the half-maximum duration of 
the precursor. In fact, the most severe restriction origi- 
nates in the condition wq (^ ~ ^b) < {'^o/b) associated 
with the dispersion correction. When 7^6, we eas- 
ily deduce from the asymptotic form of the Airy func- 
tion that the half-maximunr of the precursor will be at- 
tained for a;o(i — is) ~ 20 (wq/^)- The precursor will 
thus be well reproduced by the expression 65(2:, i) w 
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A. Square pulse 

We consider first a square-modulated incident field 
[u/f (t) — ujy(t — T)]sin(wci)- Of particular interest is 
the case where the square duration is an integer n of 
half-periods of the carrier, that is T = nTc/2 = mr/u)c- 
The incident field can then be rewritten as e(0,i) = 
unit) sin (uJct) - (-1)" unit - T) sin [w^ {t - T)] and the 
transmitted field reads as e'(z, t) = e(z, t) — (—1)" e(z, t — 
T) where e(z,t) designates the transmitted field when 
only the incident field wh (t) sin (ujct) is on. This equa- 
24000 24500 25000 tiou appllcs to thc wholc field and in particular to the 

Retarded Time (in I/Oq units) BriUouin precursor to yield: 

Figure 10: BriUouin precursor in the dispersion dominant e^(z, i) = 63(2, f) — (— 1) egyz^t — T), (37) 

limit. The solid line (bullets •) is the exact numerical so- 

lution (the analytical solution). Parameters (in cjq units): where 63(2, t) is given by Eq.lUni) or Eq.(1211), depending 
u)c = 0.836, ojp = 1.11, 7 = 10~* and ^ = 2.95 x 10"', leading on the system and the parameters considered. The two 
to cjotB ~ 2.3641 X 10", 6 ^ 0.0257cjo and 13 ^ 0.252cjo. The components of e^ are of opposite (same) sign when n is 
carrier frequency ujc is at the lower boundary of the opacity even (odd) and are well separated when it is large enough, 
region [a{uc)z si 20 J. gQ ^j^g^^ j^ significantly exceeds the duration of the ele- 

mentary precursor. On the other hand, eB{z,t) evolving 
slowly at the scale of Tc, the two components overlap and 
interfere if n is small. When n — 2 {T ~ Tc) as considered 
in [40, l45| , the two components interfere nearly destruc- 
tude if 7 < 6 and if (cjo/6)^/^ > 20, that is if b/ujQ < ^vely to give a precursor e^(z,i) « r,es(z,i - Tc./2) 



(b/ajc) Ai[—b (t ~ tg)] beyond its half-maximum ampli- 



0.024. The latter condition is approximately met FigfTUl 
for which b/ujo = 0.026. As expected, the maximum am- 
plitude of the precursor is as ~ 0.536 {h/iOc) ~ 0.0165 , 
with exp (—27 (t — tB) /3) « 0.997 at the corresponding 
time. 



VI. PROPAGATION OF PULSES WITH A 
SQUARE OR GAUSSIAN ENVELOPE 



Up to now, in the spirit of the pioneering work of Som- 
merfeld and BriUouin, we have considered incident fields 
of infinite duration. In actual or even numerical exper- 
iments, this duration is naturally finite. As a matter of 
fact the simulations made to corroborate our previous an- 
alytical calculations were made by using a square-wave 
modulation (eventually suitably filtered) and choosing a 
square duration long enough to avoid that the precursors 
generated by the rise and the fall of the square over- 
lap. On the contrary, we consider in this section the case 
where the duration of the incident field is small compared 
to the time-delay tg separating the BriUouin precursor 
from the Sommerfeld precursor and does not exceed few 
periods of the carrier. We will restrict the analysis to the 
BriUouin precursor. Indeed the Sommerfeld precursor, 
if it exists, is generally much smaller and will be often 
filtered out by rise-time effects, to which the BriUouin 
precursor is much less sensitive. 



The case where n is odd and, in particular, where n — 1 
(T — Tc/2) is much more favorable. Indeed the two pre- 
cursors then interfere constructively to yield a precursor 
e'g{z, t) w 2eB{z, t — Tc/4) whose amplitude is twice that 
obtained with a step modulation. This result is not re- 
ally a surprise since the pulse area is itself twice that of 
Unit) sin (ujct) . On the contrary the pulse area equals 
zero when n is even. The previous results are illustrated 
Figim that shows the BriUouin precursors obtained for 
n = 1 , 2 for a Lorentz medium when attenuation and dis- 
persion comparably contribute to the formation of the 
BriUouin precursor (simple asymptotic limit). 

When the detection of the BriUouin precursor is not 
time-resolved an important parameter is the integrated 



"energy" Wb{z) = j:^ \e's{z,t)\Ut 
the Parseval-Plancherel theorem [30 



[18|,|21|. Thanks to 
, it can be written 



Wb{z) 



1 

2^ 



-t-00 



expression 

|2 



\HB{z,u:)\^\E{0,u:)fduj. (38) 



all phases are eliminated and 



In this 

\Hb{z,uj)\ is reduced to exp (— 47a;^/36'^) = 
exp (— a;^/2/3^) in both strict and simple asymp- 
totic cases. For T = Tc/2, \Hb{z,uj)E{0,uj)\^ « 
(4/0;^) exp (— a;^/2/3^) and we get an energy 
WBiz) — 2'^/^7r~^/^/?/a;c which slowly decays with the 
propagation distance (as l/-y/z). On the other hand, for 
T = Tc, \Hb{z,cj)E{0,cj)\'' « {2TTUj/ojlfexp{-ojy2p^) 
and Wb{z) = (2nf/'^l3^/u}^. As expected, Wb(z) then 
decays very rapidly with the propagation distance (as 
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Figure 11: Comparison of the Brillouin precursors e'g{z,t) 
generated by an incident square-modulated field of duration 
(a) T = Tc/2 and (b) T ^ Tc . The parameters are those 
of Fig[8l The solid and dashed lines are the exact numeri- 
cal solutions, indiscernible from the analytical solutions given 
by Ea. (|37p . The bullets are the approximate solutions (a) 
2eB{z,t- Tc/l) and (b) TceB{z,t - Tc/2 . As expected the 
precursor amplitude for T = Tc/2 is twice that attained with 
a step-modulated field (see Fig[8| whereas that attained for 
T = Tc is much smaller. Inset: corresponding incident fields. 



z""^"). As already mentioned, the previous expressions 
of the energy are valid regardless of the relative con- 
tributions of the absorption and the dispersion to the 
formation of the precursor. For the Debye medium and 
the Lorentz medium in the dominant- attenuation limit, 
it is besides possible to derive from Eq. (|37|) and Eq. ([26|) 
explicit expressions of the maximum amplitude a'^ (z) of 
the precursor e'^{z, t). We find that this amplitude, equal 
to 2/3/ (wc0F) « 1.1 {P/iOc) oc 1/^2 when T = Tc/2, 
falls down to 2-\/27r/e {j3/ujc) ~ 3.0 {(3/ujc) oc 1/z when 



B. Gaussian pulse 

The Gaussian pulses are probably the sole smooth 
pulses for which it is possible to obtain exact analytic 
expressions of the Brillouin precursor, both in the strict 
and simple asymptotic limit. Non-chirped incident fields 
of the form e~* ^"^ cos (oJct) and e~* /"^ sin (wci) have 
been respectively considered by Oughstun and Balict- 
sis in [401 and by Ni and Alfano in [43|- When the 
pulses are linearly chirped, it is convenient to con- 
sider them as the real and imaginary part of e(0, t) = 
exp (iujct — t^ /T'^ + ix^t"^) where x^ is the chirping pa- 
rameter. The Fourier transform of e(0, t) and of the 
corresponding transmitted field eB{z,t) simply read as 



^(0,w) =ry^exp -{lu-ujc) T^/A 



and 



E, 



i{z,uj) = AHb{z,uj) exp (--uj^T^ /A + ujujcT^ /2] . 



(39) 



In these expressions T = T / \J\ — iy^T^ and A = 

T-yTrexp ( — a;^r^/4 j may be respectively seen as the 

(complex) duration and area of the pulse e(0,i). In the 
strict asymptotic limit [see Eq. ([M)) ]. we get 



Eb{z,uj) = ^exp 



W2 / 1 
'4! I fl2 



/32 



T2 



.T2 



(40) 



and 65(2, t), inverse Fourier transform of Eb{z,lo), reads 
as 



eB{z,t) 



Ap 



■k(i + P^t^ 



X exp 



/32 (i' _ ic.,TV2) 

1 + /32r2 



(41) 



where t' = t — ts- In the simple asymptotic limit (see 
Sec. El, Eq.dSJ) and Eq.dS!]) yield 



Eb{z,uj) — ^exp 



X exp 



—iuj \ tB + 



27 T^ 



iujrT^ 



3&3 



3ry6^ 



(42) 



This equation is easily transformed in an equation similar 
to Eq. ([32|) . By this way, we find 



(43) 



eBiz,t) =ABAi{-T]^/^bt)cxp (--rj^t 



where 



7 + 363T2 



V 



l/3fo 



exp 



-(r;V3) (27/36)^ 



and t 



B 

t - tB 



Arj'y'^/Qb'^ — iojcT"^ /2. Finally the precursors gener- 
ated by the incident fields e~* /"^ cos {ujct + x^i^) 

and e~* '^ sin {tOct + "X^t^) respectively read as 
ecos{z,t) = ^[eBiz,t)] and esin{z,t) = 3 [65(2, t)]. 
Eq. ([^T|) . Eg. (14311 and the derived expressions of ecos{z,t) 
and esin{z, t) hold whatever the duration of the incident 
pulse may be. However, as shown below, the amplitude 
of the Brillouin precursor will be only significant when 
this duration does not exceed a few periods of the 
carrier. In the Fourier transform Hb{z,uj)E{0,u>) of 
the transmitted field, Hb{z,uj) is then again much 
narrower than E{0,ll)), which may be approximated 
by its first order expansion in powers of uj. We get so 

Eb{z,uj) ^a(i + LouJcT'^/i] Hb{z,uj) and finally 



(44) 



eBiz,t) w A [hB{z,t) - (iuJcT^/2J hB{z,t) 



When there is no chirping, T and A are real, with T = T 
and A = A = TV^exp [-ujIT'^/A]. Eq.dU) then leads 
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to 



ecos{.z,t) sa AhB{z,t) 

— T^/nexp 



hB{z,t) (45) 



estn{z, t) « ^ — hsiz, t) = "—— ecos{z, t). (46) 



0.04 




Retarded Time (in 1/co„ units) 



Figure 12: Brillouin precursors generated by the incident 
fields of Gaussian envelope (a) e""'"^' cos(a;ct) with T — 

y2/cjc and (b) e"('/^'%in(a;ct) with T = VG/uJc- The pa- 
rameters are those of FigO In both cases, the pulse duration 
has been chosen in order to maximize the precursor ampli- 
tude (see text). The solid and dashed lines are the exact 
numerical solutions whereas the bullets are the analytical so- 
lutions obtained in the short pulse approximation [Ea. (|45p 
and Ea. (|46|l ]. indiscernible from those obtained without ap- 
proximation [Eq. (|43I| ]. Inset: corresponding incident fields. 
Numerical calculations shows that the Sommerfeld precursors 
generated by these fields have negligible amplitudes, respec- 
tively (a) 5.6 X 10"^ and (b) 1.15 x 10"^". 

As illustrated FigfT2l obtained in the simple asymp- 
totic limit, these approximate analytic solutions per- 
fectly fit the exact numerical solution. It is eas- 
ily deduced from Eg. pSt [Eq. p6|) ] that the ampli- 
tude of the precursor Ccosiz^t) [esi„(z,i)] is maximum 
for a pulse duration T — T„j — \/2/ujc [ VQ/(j-'c]- 
The energy of the precursors can be obtained by the 
method already used in the case of a square modu- 
lation. We get so Wb ~ (7r/2)^/^ f/ST^e-'^c^'/^^ oc 

1/yi for e(0,t) = e"*'/"^' cos (w^i) and Wb ~ 

(V32)'/' (/33w2y6g-^,^TV2^ ^ ^-3/2 fo^ g(o^^) ^ 

e~* /"^ sm{ujct). In fact the scaling laws in z^^^^ or 
z~^/^ are general and hold for every short incident pulse. 
In all cases, the transmitted pulse is indeed proportional 
to hB{z,t) when £^(0,0) = ^ ^ or to hB{z,t) when 
^ = 0, the proportionality coefficient depending only 
on the characteristics of the incident pulse and not on 
the propagation distance. For Gaussian incident pulses 



and, more generally, for smooth pulses, the amplitude 
and the energy of the Brillouin precursor rapidly de- 
creases with the pulse duration. For example, the am- 
plitude of the Brillouin precursor generated by the in- 
cident field e~* ^"^ cos (ujct) is reduced by a factor ex- 
ceeding 400 when T is taken four times larger than its 
optimum value \/2/u}c [see Eq. (|45|) ]. This reduction of 
amplitude can however be compensated by using chirped 
pulses. When the pulse duration remains small enough, 
Eq.lim holds and the Brillouin precursor generated by 
the incident field e~* /"^ cos (wc^ + X^T^) reads as 



65(2, i) 



^{z,t)^(^A^ 



hBiz,t)^(iujcAT^/2y (47) 



Anticipating that the second term of this equation is 
small compared to the first one, we easily get the ap- 
proximate expression 



(48) 



eB{z,t)KAhB z,r ~^[iujcAT^/2A 



where A = 3? 1 .4] is the area of the incident pulse. 
This result differs from that obtained without chirping 
[see Eq.(l45])[ by a extra time-delay 3? (iujcAT^/2A] and, 
moreover, by the pulse area A that may be consider- 
ably larger than that attained when the pulse is not 
chirped. Fig[T3] shows the result obtained for a pulse 
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Figure 13: Brillouin precursor generated by a chirped inci- 
dent pulse e"'-"^'^^^ cos [ujct + x^t'^) , with T = 4^/ajc and 
X ~ i^c/4 . The other parameters are as those of Fig[8] and 
Fig |12l The solid line, the bullets and the dashed line respec- 
tively are the exact numerical solution, the analytic solution 
derived from Eo. dlSp and the approximate analytic solution 
of Ea. (|48p . obtained in the short pulse approximation. Inset: 
incident pulse. The corresponding Sommerfeld precursor has 
fully negligible amplitude (9 x 10~^^ !). 

duration T ~ 4v2/wc- In order to maximize the precur- 
sor amplitude, we have chosen for the chirping the value 
X = uJc/4: for which the function A{x) reaches its first 
extremum (negative minimum). For these parameters. 



14 



5R (iAT^/2A] is also negative (time advancement). We 

remark that, despite the numerous approximations hav- 
ing led to Eq. (|48p , it provides a very good approximation 
of the exact result. 



VII. CONCLUSION 

We have analytically studied the propagation of light 
pulses in a dense Lorentz medium at distances z so large 
that the medium is opaque in a broad spectral region and 
the Sommerfeld and Brillouin precursors are far apart 
from each other. 

Assuming that the carrier frequency ujc lies in the opac- 
ity region (below, inside or beyond the anomalous disper- 
sion region) , we have shown that the Sommerfeld precur- 
sor has a shape independent of ujc and that it is entirely 
determined by the order p and the importance dp of the 
initial discontinuity of the incident field, regardless of its 
subsequent evolution. When the incident field is discon- 
tinuous (p — 0), its amplitude is independent of z and lOc ■ 
For p > 0, this amplitude is proportional to Wc^ '■^^^^■'^'' 
and rapidly decreases with the rise time of the incident 
field. These results, exact in the strict asymptotic limit 
where z — > oo, provide excellent approximations for the 
propagation distance considered by Brillouin and remain 
good approximations even when z is 1 000 times shorter. 

In the strict asymptotic limit, the formation of the Bril- 
louin precursor is uniquely determined by the frequency 
dependence of the medium attenuation. When Uc lies 
in the opacity region, we have shown that the Brillouin 
precursor is a Gaussian of amplitude as oc 1/ [ujc\/z) or 
a Gaussian-derivative of amplitude ub oc l/fw^z), de- 
pending whether the area of the incident field differs or 
not from zero. We have also determined the transmitted 
field when ujc is outside the opacity region, evidencing the 
"pollution" of the Brillouin precursor by the field that is 
then transmitted at ujc (Fig. 7). 



In a simple asymptotic limit, both attenuation and 
group delay dispersion contribute to the formation of the 
Brillouin precursor. We have established in this case an 
expression of the Brillouin precursor containing as partic- 
ular cases the previous one (dominant-attenuation limit) 
and that obtained by Brillouin by means of the stationary 
phase method (dominant-dispersion limit). 

We have finally obtained exact analytical expressions 
of the Brillouin precursors originated by pulses of square 
or Gaussian envelope. We have in particular determined 
the pulse parameters optimizing the precursor amplitude 
and demonstrated that the energy of the precursor de- 
creases with the propagation distance as slowly as z^^/^ 
when the area of the incident field differs from zero but 
as rapidly as z~^/^ in the contrary case. We have also 
shown that, for a given duration, the precursor ampli- 
tude can be greatly enhanced by using frequency-chirped 
pulses. 

Our explicit analytic expressions of the precursors con- 
trast by their simplicity from those currently derived by 
the uniform saddle point methods. The complexity of 
the latter [13j is often such that it is difficult and some- 
times impossible to retrieve from them our asymptotic 
forms. On the other hand, it should be kept in mind 
that our results only hold in the limit where the medium 
is opaque in a spectral region whose width is much larger 
than the resonance frequency. We however remark that 
they provide a not too bad reproduction of the Sommer- 
feld and Brillouin precursors even when this width is of 
the order of the resonance frequency (FiglQ]). We finally 
mention that the study of the precursors is greatly sim- 
plified when the complex index of the medium is such 
that \n{Ljj) — 1| <C 1 Vcti lla| . As in the study of the quasi- 
resonant precursors |48| . the equation giving the saddle 
points can then be reduced to a biquadratic form and the 
saddle point method is expected to provide simple solu- 
tions even when the Sommerfeld and Brillouin precursors 
partially overlap. This work is in progress. 
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